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g . ABSTRACT 

(N 

flj ' Using observations made with the XMM-Newton Observatory, we report the proba- 

Q ■ ble X-ray detection of the high- magnetic- field radio pulsar PSR J 1734-3333. This pulsar 

iy-\ . has an inferred surface dipole magnetic field of B = 5.2 x 10^'^ G, just below that of one 

anomalous X-ray pulsar (AXP). We find that the pulsar has an absorbed 0.5-2.0 keV 

flux of (5-15) X 10~^^ ergs~^ cm~^ and that its X-ray luminosity Lx is well below its 

spin down luminosity E, with Lx < 0.1^. No pulsations were detected in these data 

although our derived upper limit is unconstraining. Like most of the other high-i? pul- 

"^ ■ sars, PSR J1734-3333 is X-ray faint with no sign of magnetar activity. We tabulate 

I ' the properties of this and all other known high-i? radio pulsars with measured thermal 

^ . X-ray luminosities or luminosity upper limits, and speculate on a possible correlation 

-4— > , 

c/3 , between Lx and B. 

C^ 

^ I Subject headings: pulsars: general — pulsars: individual (PSR J1734-3333) — stars: 

^ ' neutron — X-rays: stars 

0^ 

T-H . 
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1. Introduction 
l> 

o 

f^ I Soft gamma repeaters (SGRs) and anomalous X-ray pulsars (AXPs) are now well accepted 

as being different though similar manifestations of "magnetars" — isolate d neutron sta.rs wh ose 






radiation is powered by their magnetic field (see lWoods &: Thompsonll2006l or iMereghettil 120081 for 



^ ' a review). These objects are characterised by long spin periods (2-12 s), spin-down rates that, 

^ . assuming conventional magnetic dipole braking, imply surface dipolar magnetic fields in the range 

~10^^-10^^ G, X-ray luminosities of 10^^-10^^ erg s~^ that are orders of magnitude greater than 
their spin-down luminosities, and great X-ray var iability, rang ing from short SGR-like bursts to 



major, slow-rise and long-lived X-ray flares (see iKaspil 120071 for a review). The bulk of SGR 
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and AXP properties a re well explained by the magnetar model (|Thonipson fc: Duncanlll995l . 119961 : 
Thompson et alJbood ). 



However, there remain some outstanding challenges to the magnetar picture. Arguably the 
most important is the connection between magnetars and high-magnetic field radio pulsars. Some 
models predict that given a large enough magnetic field (presumably close to t he quantum critical 
field, -Bqed = 4.4 x 10^^ G), conventional radio emission should be suppressed ([Baring fc: Harding 
19981 ) . Yet there are now known over a half dozen otherwise ordinary radio pulsars having spin-down 
inferred magnetic field B > 4 x 10^^ G (where B = 3.2 x 10^^ (PP) ' G for a pulsar with period 
P and spin-down rate P), with none showing evidence for conventional magnetar-like emission. 
Some, such as PSRs J1847-0130 and J1718-3718, have B greater than those measured for bona 
fide magnetars, yet no anomalous X-ray emission ( McLaughlin et al.ll2003l : iKaspi &: McL aughlinI 
20051) . Uncertai nties in inferred B from spin-down can be substantial (e.g., iHarding et al... 19991 : 
Spitkovskvl I2OO6I ) . but still, in the magnetar picture, some evidence for anomalous X-ray emission 
is reasonably expected in some high-i? radio pulsars. This is important overall for unifyin g the 
surprisingly varied manifestations of neutron stars into a physical framework (see lKaspill2010l for a 



review . 



Recently, the idea that high-S pulsars might exhibit anomalous X-ray emission was confirmed 
by the discovery of SGR-like X-ray bursts and a long-lived X-r ay flux enhan c emen t from what 
was previously thought to be a purely rotation-powered pulsar. iGavriil et al.l (120081 ) found that 
the high-i? pulsar PSR J1846-0258 at the center of the supernova remnant Kes 75 emitted several 
SGR-like bursts in 2006, contemporaneous with a flux enhanceme nt and a rotational glitch (see also 
Kumar &: Safi-Harbll2008l : iNg et al.ll2008l : [Livingstone et al.ll2010l ). This is the first pulsar known to 
have quiescent X-ray luminosity that could be rotation-powered, and indeed has many properties 
of rotation-powered pulsars, while showing obvious magnetar-like behavior. Its small characteristic 
age of 884 yr lends further credence to the idea — PSR J1846-0258 could be a very young magnetar, 
and one of the "missing links" in the hypothesized high-i? pulsar/magnetar evolutionary chain. 

PSR J1734-3333 is a radio pulsar with P = 1.169 s and P = 2.3 x 10~^^ which imply a 
spin-down luminosity of E = Att^IP/P^ = 5.6 x lO'^^ergs"^ and a characteristic age of r = 
P/2P = 8.1 kyr. Its inferred surface dipolar magnetic field is B = 5.2 x 10^^ G, which is among 
the highest of all known radio pulsars and just below those of bona fide magnetars such as IE 
2259+586 [B = 5.9 x 10^'^ G). From 12 yr of phase-coherent radio timing observa tions at Jodrell 
Bank , this pulsar is known to have a stable braking index measured of n = l.OibO.3 (JEspinoza et al. 



2011 



). At face value, this could imply that this high-i? pulsar's magnetic field is growing, i.e., its 
trajectory on a conventional P/P diagram is up and to the right, toward the region occupied by the 
magnetars. This makes PSR J1734-3333 a good candidate for exhibiting magnetar-like anomalous 
X-ray emission. For this reason, we obtained XMM-Newton observations of this source which we 
report on here. 



2. Observations and Results 

One 10-ks observation of PSR J1734~^333 3 was carried out on 2009 M arch 9-10 using tlie 



XMM-Newton Observatory (jjansen et alj|200ll ). The two EPIC MOS cameras (JTurner et al. 



2001 



were operating in f ull-window mode with a time resolution of 2.7 s, and the EPIC pn camera 



(jStriider et alJl200ll ) was in large-window mode providing a time resolution of 48 ms. For all three 
cameras the medium filter was in use. The total exposure time in this observation was 10.6 ks with 
each MOS camera and 8.7 ks with the pn. 

The data were analyzed with the XMM Science Analysis System (SAS) version 8.0. Co with 
calibrations updated 2009 April 13. To identify times contaminated by strong background flaring 
that is sometimes present in XMM data, we extracted and examined light curves of photons above 
10 keV over the entire field of view for the pn and MOS cameras. No flaring or significant background 
changes were found, so the entire exposure could be used in our analysis. 



2.1. Imaging 

In order to find a possible X-ray counterpart of PSR J1734-3333, we performed a blind search 
for point sources using the SAS tool edetect_chain. 

A faint X-ray source was detected near the radio position of the pulsar by edetect_chain 
in all three cameras. In the pn camera it had 49 it 10 counts in the 0.5-3.0 keV energy range 
and a likelihood ratio of L2 = — In P ~ 29 (where P is the probability that a random Poissonian 
fluctuation could produce the observed source counts). The source was also detected, in the same 
0.5-3.0 keV energy band, in the MOS 1 camera with 21 it 6 counts and L2 — 12 and in the MOS 
2 camera with 31 it 8 counts and L2 — 17. Outside this energy range there were almost no source 
counts detected. 

Figure [1] shows the X-ray emission near the radio position of the pulsar, made by combining 
the pn, MOS 1, and MOS 2 images into a mosaic and smoothing with a Gaussian kernel of radius 
a = 3" . The best-flt position of the X-ray source as reported by edetect_chain is (J2000) R.A. = 
17^34"^27!19 ± 0524, DECL. = -33°33'22'.'0 ± 3'.'0, where these uncertainties consist of both the 
1" statistical error and XMM's absolute pointmg uncertaint^H of 2". In principle it is possible to 
reduce the pointing uncertainty by matching at least two bright X-ray sources in the field to known 
optical counterparts; unfortunately, the field of view of this observation contained only one such 
source. 

The most up to date radio timing position for PSR J1734-3333 is R.A. = 17''34™26!9 ± K2, 
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DECL. = -33°33'20" ± 10" (lEspinoza et alJboilh . The offset between the radio timing position 



and the centroid of the X-ray position in R.A. is 0^29 or '^0.9 ex, and is 2'.'0 or ~0.2o" in dechnation. 



2.2. Spectroscopy 

To obtain the X-ray spectrum, we extracted a total of 300 photon events from the pn data using 
a circular region of 30" radius centered on the source. The background spectrum was taken from 
an elliptical annulus with semi- major and semi-minor axes of 75" and 45" rotated 90° from north, 
surrounding but excluding the 30" source region. This background region contained 643 events, 
implying there to be ~234 background photons in the source region. A response and ancillary 
response file were generated using the SAS tasks rmf gen and arf gen. The spectrum was grouped 
to have a minimum of 20 counts per bin with the f tool grppha and was fed into XSPEC 12.5.0ac 
for spectral fitting. We did not extract spectra from the MOS images due to the low number of 
source counts in those cameras. 

The X-ray spectrum was fit to absorbed power-law and absorbed blackbody models. Because 
of the small number of source counts, it was impossible to constrain the column density A^h- The 
Leiden/Argentine/Bonn and Dickey & Lockman Surveys of Galactic H u give the total column 
density along the line of sight to the pulsar as Ah = (1.1-1.4) x 10^^ cm^^, and from its dispersion 
measure of 578 pc cm~ ^ the distance to the pulsar is estimated to be 6.1 kpc with a ~25% uncertainty 
(jCordes &: Lazidl200ll ). We therefore assumed an upper limit on Ah of 1.2 x lO^^cm"^, halved it 
to get a more reasonable estimate of 0.6 x 10^^ cm~^ since, based on its distance, the pulsar is less 
than halfway through the Galaxy, and found best-fit models after fixing A^h to both values. The 
parameters of these best-fit models are shown in Tabled! and two of them are plotted in Figure [2j 
In general, a lower value of Ah corresponds to a power-law model with smaller photon index and 
a blackbody model with a higher kT. 



2.3. Timing Analysis 

Radio pulsar timing data were obtained using the Jodrell Bank 76-m Lovell telescope operating 
at a central obser ving frequency of aro und 1400 MHz, roughly contemporaneously with the X-ray 
observations (see lEspinoza et al.l 1201 ll for details). The TEMPO timing packagqj was used to 



co rrect the times of arriy al to the barycenter of the solar system assuming the position determined 



by lEspinoza et al.l (|201ll ) and to fit a rotational model to these arrival times. 



Based on this ephemeris for PSR J1734-3333, during the XMM observation the pulsar had 
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an expected barycentric frequency of 854.9300124 it 0.0000004 mHz. Since the MOS cameras were 
operated in full-window mode with 2.7-s time resolution, their data cannot be used for our timing 
analysis. A search for X-ray pulsations was done using the pn data after barycentering them using 
the SAS tool barycen. To do so, we extracted counts in the 0.5-3.0 keV energy range from a region 
of 20" radius centered on the source. Although we used a smaller source region here than for the 
spectral analysis in order to maximize the signal-to-noise ratio, an identical background region was 
used for both the timing and spectral work. This background region contained 286 photon events 
in the same 0.5-3.0 keV energy band, implying that, of the 84 total events in the source region, 
~46 were background photons. The source region event list was folded at the radio frequency using 
16 phase bins, and the folded light curve was fit to a constant line. The best-fit x^ was 12.2 for 15 
degrees of freedom, corresponding to a 67% probability that the folded curve could be produced 
from a data set containing no signal. Thus, no significant pulsations were detected. Searches were 
also conducted in the 0.5-10 keV and 1-2 keV energy bands with similar null results. 

To find an upper limit for the pulsed fraction, we simulated event lists with the same number 
of total counts as found in the source region. The simulated signal had a sinusoidal profile with 
a random phase and had a user specified area pulsed fraction, where the area pulsed fraction is 
defined as the ratio of the pulsed part of the profile to the entire profile. However, given the 
expected number of background counts, even a signal with 100% area pulsed fraction could only 
be detected with >3cr significance 5% of the time. Thus, we conclude that there are too few counts 
in our data set to set a meaningful upper limit on the pulsed fraction. 



3. Discussion 

Although the location of the detected X-ray source is consistent with the radio timing position 
of PSR J1734-3333, given that no X-ray pulsations were detected from it at the radio period to 
provide unambiguous proof of association, it is reasonable to question whether the X-ray source 
really is associated with the radio pulsar. We can estimate the probability of a chance superposition 
from the log A^- log S curves for XMM-Newton Galactic-plane sources in the 0.5-2.0 keV band 



(jMotchl |2006| ) . Given the lowest reasonable value of absorbed flux in that band for the source, 
~5 X 10^^^ ergs^^ cm^-^, the log A^- log S* curves predict ~180 sources per square degree at this flux 
or higher. The probability of a random X-ray source lying within the radio error ellipse is then 
only 0.1%-0.3%. 

On the other hand, looking at an optical image of the field of the XMM observation, there 
are several sources near the X-ray source. In fact, one of these optical sources (NOMAD Catalog 
ID 0564-0621454) lies within the X-ray error circle at coordinates R.A. = 17'^34™27;319 ± a;008, 
DECL. = -33°33'22'.'61 ± 0'.'05 (See Figured]). Therefore the possibility must be considered that 
the X-ray source is associated with this optical source in addition to or instead of the radio pulsar. 

Assuming that the X-ray source is associated with 0564-0621454, its X-ray to optical flux ratio 



can be estimated using the following formula from iGeorgakakis et al.1 (|2004l ): 



log(/x//< 



optj 



log / (0.5-8 keV) + 0.45 + 4.89. 



Here / (0.5-8 keV) is the unabsorbed X-ray flux in that energy band and B = 18.79 is the optical 
magnitude of 0564-0621454 taken from the USNO-Bl.O catalog. Although the unabsorbed X-ray 
flux is even less constrained than the absorbed flux used above, we can calculate a rough upper 
limit of log (/x//opt) < 1. PSR J1734--3333 is an isolated n eutron star, however, which typically 
has log(/x//opt) ~ 5 ([Treves et al.l 120001 : iKaspi et al.l 120061 ). Therefore, we conclude that if our 



X-ray source is associated with 0564-0621454, it cannot also be associated with PSR J1734-3333, 
and vice versa. 

Finally, just as the probability of a chance superposition of an X-ray source with the pulsar's 
radio position was estimated above, we can calculate the probability of an optical source coinciding 
with the X-ray position. A query of the USNO-Bl.O catalog returns almost 5500 sources within a 
10' radius of the centre of the XMM field. The probability of one of these sources lying within the 
3" X-ray error circle by chance is ~12%, much higher than the probability of a chance X-ray/radio 
alignment. Even accounting for the fact that 0564-0621454 is only 2" away from the X-ray position, 
the probability of a chance superposition is still 5%-6%. We conclude that it is highly likely that 
the X-ray source is associated with PSR J1734-3333, and much less probable that it is associated 
with the optical source 0564-0621454. 

With so few source counts to work with, the X-ray source's spectrum allows for such a wide 
range of possible photon indices and blackbody temperatures that it is impossible to eliminate a 
thermal or non-thermal source for the X-ray emission. Among the four models shown in Table [H 
the absorbed 0.5-2.0 keV flux is relatively well constrained, to within a factor of 2-3. However, 
since assuming a higher value of A'h yields a steeper power-law index or a lower kT, the unabsorbed 
0.5-2.0 keV flux is much less constrained and varies by nearly 2 orders of magnitude within and 
between the models. Assuming a distance to the pulsar of 6.1 kpc as estimated from the dispersion 
measure, if the X-ray source is the pulsar, the 0.5-10.0 keV X-ray luminosity of PSR J1734-3333 
is in the range (1-34) x lO^^ergs"^. However, Lx < E for all these models, even allowing for a 
25% uncertainty in the estimated distance, with Lx/E ~ 0.1%-10%. Furthermore, the high end of 
the luminosity range corresponds to the high A^h models, suggesting that the pulsar's true X-ray 
luminosity more likely lies in the lower end of the range, if the association holds. If it does not 
hold, the pulsar is even less luminous. 

The best-fit blackbody tem perature, kT = 0. 25 J]n'na keV, found fo r this pulsar is also unusually 



high for its age (0.07-0.11 keV; iPage et al.l 120061 ). iPons et al.l (|2007l ) reported a possible correla- 
tion between blackbody temperature and magnetic field, which, given the high S-field of PSR 
J1734-3333, could explain the discrepancy. Nevertheless, the uncertainties on kT are very large. 
Additionally we have not taken into account the effect of the neutron star atmosphere, and spectral 
models that do include the atmosphere result in a lower effective temperature. As such, without 
better statistics, we cannot eliminate the possibility that the pulsar's blackbody temperature is 



within the expected hniits. 

In Table [2] we hst all high-magnetic field radio pulsars with inferred magnetic field B > 
-SqeDi as well as all other radio pulsars with B > 1.5 x 10^'^ G, that have measured thermal X-ray 
luminosities or luminosity up per limits. Of thes e stars only one, PSR J1846-0258, has exhibited 



clear magnetar-like behavior (iGavriil et al.ll2008l ). Other than that, the high-i? pulsars show X-ray 
spectral properties that are much diff'erent from those of active magnetars, being much softer and 
fainter X-ray sources. They do, however, show some similarities with transient magnetars, which 
in quiescence are also softer, fainter X-ray sources. One transient magnetar in particular, XTE 
J1810-197, displays spectral properties in quiescence that are consistent with those of the de tected 



high-S pulsars {kT « 0.18 keV and 0.5-10 keV Lx ~ 7 x 10^^ ergs'^) (jGotthelf et al.l 12004 1 



Additionally, we took the reported X-ray luminosities for the pulsars in Table El converted them 
to the 0.5-10 keV energy band, and plotted this value versus inferred magnetic field in Figure [3j 
The resulting plot shows possible evidence of a correlation between the two quantities, with higher 
S-field stars having higher 0.5-10 keV Lx', certainly the pulsars here with i? > 4 x 10^^ G are 
brighter in X-rays than those with i? < 4 x 10^^ G. There are important caveats that should be 
noted, however. In particular, the four youngest pulsars in Table [2] also comprise four out of the 
five brightest X-ray detected pulsars in the Table, suggesting that the perceived correlation may 
be due entirely or in part to age instead of B. On the other hand, one of the brightest of these 
pulsars, PSR J1819~1458, has a considerably higher characteristic age; it is higher even than that 
of the much dimmer PSR B1916-I-14. Additionally, of the four pulsars with the highest B fields, 
two have only X-ray luminosity upper limits and the other two have only weak detections with 
poorly constrained luminosity measurements. Future observations may show some or all of these 
pulsars to be substantially more X-ray dim than PSR J1119-6127 or PSR J1819~1458, which would 
considerably weaken the possible correlation. 

Although the spectrum of PSR J1734~3333 is poorly constrained, it is evident that it has much 
more in common with the other high-S rotation-powered pulsars than with the magnetars. Most 
magnetars, even in quiescence, have X-ray luminosities exceeding their spin-down luminosities 
by orders of magnitude; for this pulsar, as with the other high-S radio pulsars, the opposite is 
true. In particular, its 2-1 keV X-ray luminosity is well below that of any known magnetar (see 



Woods fc Thompson! I2OO6I . Table 14.1). Therefore, whil e PSR J1734-3333 co uld be argued to be 



evolving towards the rotational properties of a magnetar (JEspinoza et al.ll201lh . it currently retains 
essentially the X-ray radiative properties of a radio pulsar. What is less clear from these data, 
however, is whether its 0.5-10 keV Lx/E is similar to or greater than that of PSR J1846-0258 in 
quiescence (~1%-10%), or whether it is considerably smaller as with most other high-B pulsars. 
Unfortunately, due to the poor statistics it is difficult to compare the spectrum of PSR J1734-3333 
with those of other high-i? pulsars and magnetars. Longer X-ray observations are needed to better 
constrain the luminosity, obtain a good spectrum, and search for X-ray pulsations. 
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Fig. 1.— XMM image of the PSR J1734-3333 field in the 0.5 - 3.0 keV band, smoothed by a 
Gaussian kernel with a = 3". The radio timing position is shown by the ellipse, and the crosses 
denote the positions of optical sources from the USNO-Bl.O catalog. The optical source closest to 
the X-ray source (NOMAD Catalogue ID 0564-0621454) is represented by the cross marked with a 
box. 
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Fig. 2. — XMM pn spectrum of the possible X-ray counterpart to PSR J1734-3333. The spectrum 
is binned to contain a minimum of 20 counts per bin. Top: best-fit blackbody model for fixed 
A^'h = 0.6 X 10^^ cm^^. Bottom: best-fit power-law model for fixed A'^h = 1-2 x 10^^ cm^^. 
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Fig. 3. — X-ray luminosity (0.5-10 keV) vs. inferred surface dipolar magnetic field strength for X- 
ray observed high-S radio pulsars. Luminosity values from the original references were extrapolated 
to this energy band based on the fit or assumed spectral model. Error bars are at the la confidence 
level and do not include uncertainties in the distance measurements. Note that the value plotted 
for PSR J1846-0258 is its quiescent X-ray luminosity. The dotted line represents the quantum 
critical field, -Bqed = 4.4 x 10^^ G. See Table [2] for references. 
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Table 1. Spectral Parameters for PSR J1734-3333 



Parameter 




Power Law 


Blackbody 


Nh (1022 cm-2) 




1.2 


0.6 


1.2 


0.6 


pa 




4.2+}i 


9 fi+i-3 






kT (keV)^ 








n 9f;+0-i3 

'-'•^^-0.08 


0.34+°:}? 


X' (dof) 




1.44(7) 


2.55(7) 


1.35(7) 


1.49(7) 


/abs (lO-i^ergs-icm- 


r 


6.4-16 


5.8-15 


5.7-15 


6.1-15 


/unabs (iQ-i'^ergs^icm 


-2y 


6.7-77 


1.7-7.5 


3.6-21 


1.5-4.9 



^Uncertainty ranges indicate 90% confidence intervals. 

'^Absorbed flux in the 0.5-2.0 keV band. Range indicates 90% confi- 
dence interval. 

'^Unabsorbed flux in the 0.5-2.0 keV band. Range indicates 90% con- 
fidence interval. 



Table 2. High-magnetic- field Radio Pulsars 



Name 


P(s) 


B (1013 G) 


E (ergs-i) 


r (kyr) 


d" (kpc) 


Lx"^ (ergs-i) 


Energy Range (keV) 


Reference 




J1847-0130 


6.71 


9.4 


1.7 X 


10^2 


83 


8.4 


< 5 X 10^3 


2-10 


McLaughlin et al. ("2003) 




J1718-3718 


3.38 


7.4 


1.6 X 


10^3 


34 


4.5 


0.2-6 X 10^3 


Bolometric 


Kaspi & McLaughlin (2005) 




J1814-1744 


3.98 


5.5 


4.7 X 


1032 


85 


10 


< 6.3 X 10^5 

< 4.3 X 10^3 


0.1-2.4 
2-10 


Pivovaroff et al. (2000) 




J1734-3333 


1.17 


5.2 


5.6 X 


10^4 


8.1 


6.1 


0.1-3.4 X 10^3 
0.03-2.2 X 10^^ 


0.5-10 
2-10 


This work 




J1819-1458'= 


4.26 


5.0 


2.9 X 


1032 


117 


3.6 


2.8-4.3 X 10^3 


0.3-5 


McLaughhn et al. (2007) 




J1846-0258 


0.33 


4.9 


8.1 X 


1036 


0.9 


6.0'' 


2.5-2.8 X IQ^" 
1.2-1.7 X lO^'^'^ 


0.5-10 
0.5-10 


Ng et al. (2008) 
Ng et al. (2008) 




J1119-6127 


0.41 


4.1 


2.3 X 


1036 


1.7 


8.4<= 


1.9-3.2 X 10^3 


0.5-7 


Safi-Harb & Kumar (2008) 




J0847-4316" 


5.98 


2.7 


2.2 X 


10^1 


790 


3.4 


< 1 X 10^2 


0.3-8 


Kaplan et al. (2009) 




J1846-0257'= 


4.48 


2.7 


7.1 X 


10^1 


442 


5.2 


< 3 X 10^2 


0.3-8 


Kaplan et al. (2009) 




B0154+61 


2.35 


2.1 


5.7 X 


1032 


197 


1.7 


< 8 X 10^1 


0.3-10 


Gonzalez et al. (2004) 


1 


B1916+14 


1.18 


1.6 


5.1 X 


1033 


88 


2.1 


~ 3 X 10^1 


Bolometric 


Zhu et al. (2009) 


(J^ 



''Unless otherwise noted, all distances were estimated from the dispersion measure of the source. 

•"The ranges for PSRs J1718-3718 and J1819-1458 indicate 68% confidence intervals, while the ranges for PSRs J1734-3333, J1846-0258, and 
J1119-6127 indicate 90% confidence intervals. 

"^Pulsar classified as a rotating radio transient (RRAT). 

"^The distance to PSR J1846-0258 was found from H I and '^CO spectral measurements JLeahv fc Tianll2008l ). 



''The distance to PSR J1119-6127 was found from H I absorption measurements IjCaswell et al.ll2004l ) 
This value is from 2000, prior to this pulsar's 2006 outburst. 
^This value is from 2006, during this pulsar's 2006 outburst. 



